Expression and high levels of insertional polymorphism of an endogenous gammaretrovirus lineage in dogs.

Despite the absence of a confirmed exogenously replicating retrovirus in Canis lupus familiaris (C. familiaris), past retroviral infections are evident in the genomes of living animals via the presence of endogenous retroviruses (ERVs). Although gammaretrovirus-like transcripts and enzyme activities were previously reported to be present in canine leukemias and lymphomas, those findings were not further explored. Initial analysis of the C. familiaris reference genome revealed a minor subset of one ERV lineage, classified as CfERV-Fc1(a), or Fc1(a) here, with features characteristic of recent integration, including the presence of ORFs and identical or nearly identical LTRs. Our previous analysis of whole genome sequence data belonging to extant Canidae revealed a burst of past infections in Canis ancestors resulting in numerous young, polymorphic, and highly intact loci now segregating in dogs. Here, we demonstrate the expression of full-length Fc1(a) proviruses in tissues collected from healthy animals and from animals with cancer. We observed significantly higher expression in samples of dogs with various cancer diagnoses when compared to samples from healthy dogs. Genotyping of insertionally polymorphic Fc1(a) loci identified candidate expressed proviruses and delineated distributions over sample groups. Collectively, the data show that Fc1(a) proviruses retain biological activity in the domestic dog and provides a means to examine potential genetic links with disease states in this species.

A reference genome for the long-term kleptoplast-retaining sea slug Elysia crispata morphotype clarki.

Several species of sacoglossan sea slugs possess the incredible ability to sequester chloroplasts from the algae they consume. These "photosynthetic animals" incorporate stolen chloroplasts, called kleptoplasts, into the epithelial cells of tubules that extend from their digestive tracts throughout their bodies. The mechanism by which these slugs maintain functioning kleptoplasts in the absence of an algal nuclear genome is unknown. Here, we report a draft genome of the sacoglossan slug Elysia crispata morphotype clarki, a morphotype native to the Florida Keys that can retain photosynthetically active kleptoplasts for several months without feeding. We used a combination of Oxford Nanopore Technologies long reads and Illumina short reads to produce a 786-Mb assembly (N50 = 0.459 Mb) containing 68,514 predicted protein-coding genes. A phylogenetic analysis found no evidence of horizontal acquisition of genes from algae. We performed gene family and gene expression analyses to identify E. crispata genes unique to kleptoplast-containing slugs that were more highly expressed in fed versus unfed developmental life stages. Consistent with analyses in other kleptoplastic slugs, our investigation suggests that genes encoding lectin carbohydrate-binding proteins and those involved in regulation of reactive oxygen species and immunity may play a role in kleptoplast retention. Lastly, we identified four polyketide synthase genes that could potentially encode proteins producing UV- and oxidation-blocking compounds in slug cell membranes. The genome of E. crispata is a quality resource that provides potential targets for functional analyses and enables further investigation into the evolution and mechanisms of kleptoplasty in animals.

Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote.

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application.

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

Reference Genome for the Highly Transformable Setaria viridis ME034V.

Setaria viridis (green foxtail) is an important model system for improving cereal crops due to its diploid genome, ease of cultivation, and use of C4 photosynthesis. The S. viridis accession ME034V is exceptionally transformable, but the lack of a sequenced genome for this accession has limited its utility. We present a 397 Mb highly contiguous de novo assembly of ME034V using ultra-long nanopore sequencing technology (read N50 = 41kb). We estimate that this genome is largely complete based on our updated k-mer based genome size estimate of 401 Mb for S. viridis Genome annotation identified 37,908 protein-coding genes and >300k repetitive elements comprising 46% of the genome. We compared the ME034V assembly with two other previously sequenced Setaria genomes as well as to a diversity panel of 235 S. viridis accessions. We found the genome assemblies to be largely syntenic, but numerous unique polymorphic structural variants were discovered. Several ME034V deletions may be associated with recent retrotransposition of copia and gypsy LTR repeat families, as evidenced by their low genotype frequencies in the sampled population. Lastly, we performed a phylogenomic analysis to identify gene families that have expanded in Setaria, including those involved in specialized metabolism and plant defense response. The high continuity of the ME034V genome assembly validates the utility of ultra-long DNA sequencing to improve genetic resources for emerging model organisms. Structural variation present in Setaria illustrates the importance of obtaining the proper genome reference for genetic experiments. Thus, we anticipate that the ME034V genome will be of significant utility for the Setaria research community.

Conservation of Cdc14 phosphatase specificity in plant fungal pathogens: implications for antifungal development.

Cdc14 protein phosphatases play an important role in plant infection by several fungal pathogens. This and other properties of Cdc14 enzymes make them an intriguing target for development of new antifungal crop treatments. Active site architecture and substrate specificity of Cdc14 from the model fungus Saccharomyces cerevisiae (ScCdc14) are well-defined and unique among characterized phosphatases. Cdc14 appears absent from some model plants. However, the extent of conservation of Cdc14 sequence, structure, and specificity in fungal plant pathogens is unknown. We addressed this by performing a comprehensive phylogenetic analysis of the Cdc14 family and comparing the conservation of active site structure and specificity among a sampling of plant pathogen Cdc14 homologs. We show that Cdc14 was lost in the common ancestor of angiosperm plants but is ubiquitous in ascomycete and basidiomycete fungi. The unique substrate specificity of ScCdc14 was invariant in homologs from eight diverse species of dikarya, suggesting it is conserved across the lineage. A synthetic substrate mimetic inhibited diverse fungal Cdc14 homologs with similar low µM Ki values, but had little effect on related phosphatases. Our results justify future exploration of Cdc14 as a broad spectrum antifungal target for plant protection.

Hybrid de novo genome assembly of red gromwell (Lithospermum erythrorhizon) reveals evolutionary insight into shikonin biosynthesis.

Lithospermum erythrorhizon (red gromwell; zicao) is a medicinal and economically valuable plant belonging to the Boraginaceae family. Roots from L. erythrorhizon have been used for centuries based on the antiviral and wound-healing properties produced from the bioactive compound shikonin and its derivatives. More recently, shikonin, its enantiomer alkannin, and several other shikonin/alkannin derivatives have collectively emerged as valuable natural colorants and as novel drug scaffolds. Despite several transcriptomes and proteomes having been generated from L. erythrorhizon, a reference genome is still unavailable. This has limited investigations into elucidating the shikonin/alkannin pathway and understanding its evolutionary and ecological significance. In this study, we obtained a de novo genome assembly for L. erythrorhizon using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies. The resulting genome is ∼367.41 Mb long, with a contig N50 size of 314.31 kb and 27,720 predicted protein-coding genes. Using the L. erythrorhizon genome, we identified several additional p-hydroxybenzoate:geranyltransferase (PGT) homologs and provide insight into their evolutionary history. Phylogenetic analysis of prenyltransferases suggests that PGTs originated in a common ancestor of modern shikonin/alkannin-producing Boraginaceous species, likely from a retrotransposition-derived duplication event of an ancestral prenyltransferase gene. Furthermore, knocking down expression of LePGT1 in L. erythrorhizon hairy root lines revealed that LePGT1 is predominantly responsible for shikonin production early in culture establishment. Taken together, the reference genome reported in this study and the provided analysis on the evolutionary origin of shikonin/alkannin biosynthesis will guide elucidation of the remainder of the pathway.

Origin and recent expansion of an endogenous gammaretroviral lineage in domestic and wild canids.

BACKGROUND: Vertebrate genomes contain a record of retroviruses that invaded the germlines of ancestral hosts and are passed to offspring as endogenous retroviruses (ERVs). ERVs can impact host function since they contain the necessary sequences for expression within the host. Dogs are an important system for the study of disease and evolution, yet no substantiated reports of infectious retroviruses in dogs exist. Here, we utilized Illumina whole genome sequence data to assess the origin and evolution of a recently active gammaretroviral lineage in domestic and wild canids. RESULTS: We identified numerous recently integrated loci of a canid-specific ERV-Fc sublineage within Canis, including 58 insertions that were absent from the reference assembly. Insertions were found throughout the dog genome including within and near gene models. By comparison of orthologous occupied sites, we characterized element prevalence across 332 genomes including all nine extant canid species, revealing evolutionary patterns of ERV-Fc segregation among species as well as subpopulations. CONCLUSIONS: Sequence analysis revealed common disruptive mutations, suggesting a predominant form of ERV-Fc spread by trans complementation of defective proviruses. ERV-Fc activity included multiple circulating variants that infected canid ancestors from the last 20 million to within 1.6 million years, with recent bursts of germline invasion in the sublineage leading to wolves and dogs.

Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication.

BACKGROUND: Domesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication. RESULTS: Whole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication. CONCLUSIONS: Our results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.

Ancient European dog genomes reveal continuity since the Early Neolithic.

Europe has played a major role in dog evolution, harbouring the oldest uncontested Palaeolithic remains and having been the centre of modern dog breed creation. Here we sequence the genomes of an Early and End Neolithic dog from Germany, including a sample associated with an early European farming community. Both dogs demonstrate continuity with each other and predominantly share ancestry with modern European dogs, contradicting a previously suggested Late Neolithic population replacement. We find no genetic evidence to support the recent hypothesis proposing dual origins of dog domestication. By calibrating the mutation rate using our oldest dog, we narrow the timing of dog domestication to 20,000-40,000 years ago. Interestingly, we do not observe the extreme copy number expansion of the AMY2B gene characteristic of modern dogs that has previously been proposed as an adaptation to a starch-rich diet driven by the widespread adoption of agriculture in the Neolithic.

Duplications and losses in gene families of rust pathogens highlight putative effectors.

Rust fungi are a group of fungal pathogens that cause some of the world's most destructive diseases of trees and crops. A shared characteristic among rust fungi is obligate biotrophy, the inability to complete a lifecycle without a host. This dependence on a host species likely affects patterns of gene expansion, contraction, and innovation within rust pathogen genomes. The establishment of disease by biotrophic pathogens is reliant upon effector proteins that are encoded in the fungal genome and secreted from the pathogen into the host's cell apoplast or within the cells. This study uses a comparative genomic approach to elucidate putative effectors and determine their evolutionary histories. We used OrthoMCL to identify nearly 20,000 gene families in proteomes of 16 diverse fungal species, which include 15 basidiomycetes and one ascomycete. We inferred patterns of duplication and loss for each gene family and identified families with distinctive patterns of expansion/contraction associated with the evolution of rust fungal genomes. To recognize potential contributors for the unique features of rust pathogens, we identified families harboring secreted proteins that: (i) arose or expanded in rust pathogens relative to other fungi, or (ii) contracted or were lost in rust fungal genomes. While the origin of rust fungi appears to be associated with considerable gene loss, there are many gene duplications associated with each sampled rust fungal genome. We also highlight two putative effector gene families that have expanded in Cqf that we hypothesize have roles in pathogenicity.